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Density functional theory calculations have been used to identify stable layered Li-M-B crystal 
structure phases derived from a recently proposed binary metal-sandwich (MS) lithium monoboride 
superconductor. We show that the MS lithium monoboride gains in stability when alloyed with 
electron-rich metal diborides; the resulting ordered Li 2 (i- X )M X ~B2 ternary phases may form under 
normal synthesis conditions in a wide concentration range of x for a number of group-III-V metals 
M. In an effort to pre-select compounds with the strongest electron-phonon coupling we examine 
the softening of the in-plane boron phonon mode at F in a large class of metal borides. Our results 
reveal interesting general trends for the frequency of the in-plane boron phonon modes as a function 
of the boron-boron bond length and the valence of the metal. One of the candidates with a promise 
to be an MgB2-type superconductor, Li2AlB4, has been examined in more detail: according to our 
ab initio calculations of the phonon dispersion and the electron-phonon coupling A, the compound 
should have a critical temperature of ~ 4 K. 



I. Introduction 



Observation and explanation of the superconducting 
transition in MgB2 at an unexpectedly high temperature 
of 39 K [l|, |2[ have stimulated extensive research aimed 
at developing related layered phonon- mediated supercon- 
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Interesting 



ductors [: 

results have been recently obtained for carbon-based lay- 
ered materials: the efforts to adjust the properties of 
graphite via intercalation with various metals have led 
to the discovery of a CaCs sup erconductor with a critical 
temperature of 11.5 K [jj llOj (15.1 K under hydrostatic 
pressure fill]). 

Tuning the properties of boron-based superconductors 
has proven to be difficult because, despite the existence 
of over a dozen of stable metal diborides with the C32 
only one of them, 



structure 



MgB 2 , has the 
hole-doped boron pa states that couple strongly to the 
in-plane boron phonon modes For comparison, crit- 
ical temperatures in other C32-MB 2 superconductors, 
such as TaB 2 or NbB 2 0, H, @> do not exceed 10 K be- 
cause these electron-rich compounds have fully occupied 
boron pa bands and thus lack the important nearly two- 
dimensional Fermi surfaces. The attempts to raise T c in 
MgB 2 via substitutional doping have been unsuccessful 
for various reasons: i) electron-doping led to the decrease 
of the electron-phonon coupling due mainly to the filling 
of the pa band Q ; ii) hole-doping or substitution of Mg 
with large metals turned out to be thermodynamically 
unfavorable 17 , li| ; and iii) co-doping of MgB 2 with Li- 
Al Li-C caused the reduction of the T c because of the 
concurrent filling of the pa band and depletion of the p7r 



band |l9j. Recently, Palnichcnko et al. have reported 
interesting results on a possible onset of superconductiv- 
ity in MgB 2 after thermal treatment in the presence of 
Rb, Cs, and Ba; the structure of the resultant materials 
and the mechanism of the T c enhancement have yet to 
be determined ^EH. 

Identification and synthesis of new stable compounds 
will be a critical step to overcome the limitations of 
the existing layered metal borides and to have a chance 
of obtaining better superconductors. Using a data- 
mining approach we have recently found a previously 
unknown metal sandwich (MS) crystal structure and 
demonstrated that lithium monoboride in this configu- 
ration is marginally stable under ambient conditions and 
is likely to form under pressure 
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Remarkably, 

our calculations indicated that, relative to MgB 2 , MS- 
LiB has a higher density of boron pa states at the Fermi 
level 
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22|, a feature long sought for in MgB 2 -related 
materials. Subsequent theoretical studies of the electron- 
phonon coupling predicted that the T c in MS-LiB would 
be around 10 K [S HI- The relatively low value of T c 
has been attributed to an accidental electronic structure 
feature present in the pristine MS-LiB: the crossing of 
the pn bands happens to be exactly at the Fermi level 



23l . |24| . After examining the compound's response to 



hydrostatic pressure and substitutional doping we have 
concluded that the recovery of the pir density of states 
(DOS) at the Fermi level would be accompanied by an 
unfortunate drop of the pa DOS under these conditions 
[24j |. The question of whether it could be possible to 
tune the properties of MS-LiB to obtain higher T c re- 
mains open. 

In this work we expand the search for stable crystal 



2 



structure phases to ternary metal borides. Namely, we 
use density functional theory calculations to explore the 
relative stability of different layered configurations with 
composition l^i2(i-x)M x B 2 for a library of over 30 met- 
als M. We show that MS lithium monoboride gains in 
stability when alloyed with electron-rich metal diborides; 
the resulting ordered Li 2 (i_ a; )M a ;B2 ternary phases may 
form under normal synthesis conditions in a large concen- 
tration range of x for a number of group-III-V metals M. 
For several promising candidates we evaluate the soften- 
ing of the in-plane boron phonon mode at Y in order to 
identify which of these compounds has large coupling to 
the in-plane boron modes (similarly to what happens in 
MgB2). For one case, L^AlB-i, we also perform a full 
calculation of the critical temperature. We believe that 
our results will be of considerable help to experimental 
groups working on the development of new boron-based 
layered superconductors. 

The paper is organized as follows. In Sec. II we de- 
scribe the geometry of the ternary MS structures. Section 
III is devoted to the zero-temperature thermodynamic 
stability of IA 2 (\- X )M X Q2 phases and to the special case 
Li2(i-x)ALt32. Electronic properties of two representa- 
tives, Li2AlB4 and Li2TiB4, are addressed in Sec. IV. 
In Sec. V we examine the phonon softening in a large 
set of metal borides. In Sec. VI we present the ab initio 
superconductivity analysis of Li2AlB4. Conclusions are 
given in Sec. VII. 
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FIG. 1: (Color online). Structure of the proposed ternary 
metal sandwich phases and their location on a schematic 
ternary (Li-M-B) phase diagram. Formation energies are 
shown along the vertical axis. The relative stability is cal- 
culated with respect to 11-LiB (MS2) and 0-MB 2 (C32) (no- 
tation explained in the text). 
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M X ~B2 metal sandwich structures 



Consideration of ternary Li-M-B systems dramatically 
increases the number of candidates to be screened be- 
cause one should check both different compositions and 
various populations of metal sites. We reduce the search 
space by recalling the binding trends for the binary MS 
metal borides demonstrated in our previous study [22| . 
Namely, the MS lithium monoboride has plenty of avail- 
able bonding boron pa states and may additionally sta- 
bilize when mixed with electron-rich metals. Because 
straightforward substitutional doping generally causes 
strain and loss of binding within the metal layers ID, 25 1 
we separate regions with different metal species by boron 
layers and, thus, avoid putting different metals in direct 
contact. This idea follows a natural segregation tendency 
in layered borides, which has been observed, for instance, 
in MgB 2 under heavy Al doping: the resulting compound 
was a C32 superstructure with alternating Mg and Al 
layers firj] ]. 

Our candidate segregated superstructures are built by 
combining MS-LLB and C32-MB2 unit cells. Since there 
are numerous possibilities to stack metal and boron lay- 



ers, it is convenient to adopt the following compact nota- 
tion: MSl-LiB and C32-AfB 2 basic units will be denoted 
as "1"- and "0"-unit, respectively. In this way one can 
uniquely specify the structure with a string of ones and 
zeros. Note that the presence of one "l"-unit results in 
a stacking shift along the c-axis which makes the unit 



cell rhombohedral 12, k22j|. For convenience, in candi- 



date phases containing an even number of "1" -units we 
choose to reflect every second "l"-unit in the (x, y)-plane 
to obtain a hexagonal unit cell. This choice is expected 
to not affect our results, since we have found no appre- 
ciable differences in the properties of the rhombohedral 
and hexagonal allotropes of MS-LiB [12]. Examples of 
two superstructures with small unit cells are 10-Li2MB4 
[Fig. W[c)] and 110-Li 4 MB 6 [Fig. [Ifb)], the latter being 
the superposition of MS2-LiB and C32-MB 2 . 

Figure Q] shows that such compounds are located along 
the red line on the ternary phase diagram. Note that 
MS-LiB is marginally stable (the known competing phase 
LiBj, has a completely different structure |22j). This im- 
plies that if (i) the C32-MB2 exists and (ii) the relative 
stability of the MS ternary phases along this line is neg- 
ative, then one would have a better chance of synthesiz- 
ing the MS ternary compounds from elemental materials 
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rather than the two binaries MS2-LiB and MB2 sepa- 
rately. Of course, one needs to make sure that no other 
competing phase forms, which can be a challenging task 
for ternary systems. 



III. Thermodynamic stability 

Computational details. In the evaluation of the low 
temperature stability of ternary compounds, we rely on 
the " ah initio formation energy" criterion, which has 
been shown to be a reliable approach for binary systems 
[27L [28j| (expected probability of predicting the correct ex- 
perimental compounds or immiscibilities is r\* « 97.3% 
as defined in equation (5) of Ref. [29[). Ab initio to- 
tal energy calculations are performed in the generalized 
gradient approximation with projector augmented-wave 
pseudopotentials (PAW) [3(| and Perdew, Burke, and 
Ernzerhof [3 11 ] exchange-correlation functional, as imple- 
mented in VASP [32I • Because of a significant charge 
transfer between metal and boron in most structures con- 
sidered we use PAW pseudopotentials in which semicore 
states are treated as valence. This is especially impor- 
tant for the Li-B system as discussed in Rcfs. [iH and 
. Simulations are carried out at zero temperature and 
without zero-point motion; spin polarization is used only 
for borides of Fe, Co, and Ni. We use an energy cutoff of 
398 eV and at least 4000/ (number of atoms in unit cell) 
k-points distributed on a Monkhorst-Pack mesh [35| . All 
structures are fully relaxed and numerically converged to 
within 1-2 meV/atom. 

Stability of Li%MB± phases. Figure [2] shows the rela- 
tive formation energy AEj° of 10-Li2-A/fB 4 metal borides 
(Fig. [He)) with respect to phase separation into 11-LiB 
(MS2) and 0-MB 2 (C32): 
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(all formation energies are per atom). We observe that 
there are many metals that can stabilize MS2-LiB and 
a number of them do have stable C32-MB2 phases (M 
= Al, Hf, Ti, V, Nb, Ta) [H, El- Considering that 
metals M — Hf, Ti, V, Nb, Ta have no reported sta- 
ble compounds with Li or Li-B [l5j . there is a good 
chance that the predicted layered phases will form in 
the Li-M-B (M = V, Nb, Ta) ternary systems. Note 
that stable structures for several transition-metal di- 
borides have larger unit cells with buckled boron lay- 
ers and different location of metal sites [WB2 (hP12, 
hR18), MoB 2 (hR18), RuB 2 and OsB 2 (0P6) (the 5 
phase [37| is basically a non-corrugated 0P6)]. For ex- 
ample, the 0P6-RUB2 phase is favored over C32-RuB 2 
by about 0.4 eV/atom, which would make 10-Li2RuB4 
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FIG. 2: (Color online). Relative stability of the 10-Li 2 A/B 4 
metal borides (Fig. QJc)) with respect to phase separation 
into 11-LiB (MS2) and O-MB2 (C32) (see Eq. 1) as a function 
of metal valence. The group-II-V metals highlighted in blue 
form stable C32-type diborides. 



unstable. However, one could expand the library of pos- 
sible ternary configurations by creating more stable MS 
sequences with buckled boron layers. As for the mono- 
and divalent metal diborides, it is not surprising that 
they cannot improve the stability of MS-LiB in the cor- 
responding ternary alloy because they are underdoped 
themselves 16]. Structural and electronic properties for 
selected proposed 1010-Li2-MB4 phases are summarized 
in Table I. 

Li-Al-B system. Aluminum is a special case because 
it is the least electron-rich metal on the list that pro- 
vides the desired stabilization for the ternary MS con- 
figurations. The sizeable 50-meV/atom energy gain 
for 10-Li2AlB4 indicates that this compound would be 
stable with respect to C32-A1B 2 and the known off- 
stoichiometry LiB y as well [22[ ■ For the analysis of the 
thermodynamic stability of Li 2 ( x _i) A1 X B2 we also need 
to consider the known stable Li-Al binary phases B32- 
LiAl, hR15-Li 3 Al 2 , and mS26-Li 9 Al 4 [ll,[Mpj| and the 
only reported ternary LiAlBi4 compound 40] which was 
observed in experiments on doping AIB2 with Li [411 ]. 

The location of these phases on the Li-Al-B phase di- 
agram is shown in Fig. 3. Because 10-Li2AlB4 lies very 
close to the line connecting LiAlBi4 and mS26-LigAl 4 , 
its relative stability will be defined primarily by the for- 
mation energies of these two compounds. We find that 
10-Li2AlB4 is virtually degenerate in energy with respect 
to the mixture of LiAlBi4, mS26-LigAl4, and Li at zero 
temperature and pressure and without zero-point correc- 
tions; the finite temperature contributions must have a 
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TABLE I: Structural and electronic properties of the proposed 1010-Li2-MB4 phases for metals M that have stable C32-MB2 
phases. We list fully relaxed lattice vectors and Wyckoff positions for the PGz/mmc (no. 194) unit cells: Bl (4/) (1/3,2/3,zbi); 
B2 (4e) (0,0, z B2 ); Li (4/) (l/3,2/3,* w ); and M (2c) (1/3,2/3,1/4). E^ 1 and E^' 2 and positions of the boron a states above 
the Fermi level at the F point and Ej 010 is a relative formation energy per atom with respect to MS2-LiB and C32-MB2. Band 
structures for 1010-Li2AlB4 and 1010-Li2TiB4 are shown in Figs. 5 and 6, respectively. 
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negative effect on the 10-Li2AlB4 relative stability, as 
there are no reports on the formation of this compound. 

In order to stabilize the proposed phase one could use 
high pressures: as we have demonstrated previously the 
MS phases are unusually soft [l2j, l22j. Indeed, the calcu- 
lated relative formation enthalpy as a function of pres- 
sure for 10-Li2AlB4 in Fig. 3 either becomes negative 
(with respect to LiAlBi4, mS26-LigAi4, and Li) or re- 
mains negative (with respect to 11-LiB and O-AIB2). The 
considerable change in slope at about 6 GPa is related 
to a sudden ~ 10% decrease in the 10-Li 2 AlB4 atomic 
volume at that pressure (Fig. 3). Overall, 10-Li2AlB4 
compresses by over 25% when the pressure is increased 
from GPa (v = 10.4 A 3 /atom, H f = -0.163 eV/atom) 
to 12 GPa (v = 7.7 A 3 /atom, H f = -0.302 eV/atom). 
For comparison, the LiAlBi4 phase is much more com- 
pact at P = GPa (v = 7.7 A 3 /atom, H f = -0.168 
eV/atom) but compresses by only 5% at P = 12 GPa 
(v = 7.3 A 3 /atom, H f = -0.217 eV/atom). These re- 
sults suggest that this ternary and other binary boron- 
rich Li-B and Al-B phases [lilli^ , also comprised of rigid 
boron polyhedra (e.g., L13B14), should not prohibit the 
formation of the predicted layered compounds under the 
pressures considered. 

There is strong evidence that another stable phase has 
formed in different experiments involving Li, Al, and B 
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43j. Monni et al. suggested that their unidentified 
low- angle spurious peaks could belong to a binary Li-B 
phase (43j. The authors purposely prepared a sample at 
the LiB2 composition and indeed observed a matching 
set of the low- angle reflections (see Fig. 7 in Ref. [43]). 
Their conclusion that the unknown phase is a binary Li-B 



compound agrees well with the observation of 12.2° and 
20.9° peaks at 40% - 50% Li compositions by Wang et al. 
in 1978 [44J. However, in a series of detailed studies of the 
Li-B system Worle et al. have conclusively determined 
that such low-angle peaks correspond to Li 6 B 18 (Li 2 0) 2 ; 
phases (0< x <1) [45j. The authors' latest analysis of the 
X-ray patterns and electron density distributions has in- 
dicated that the zeolite-like structure of LigBi8(Li 2 0) x 
consists of interconnected boron octahedra with large 
tunnels filled with neutral Li2 template [13] . There- 
fore, it appears that the stable oxygen-containing phase 
also formed in the earlier experiments involving Li and B 
at about 1:2 composition. This knowledge is important 
for future experiments on the ternary metal borides: in 
order to obtain the predicted layered compounds (whose 
stability can be affected by a few meV/atom) one needs 
to reduce the amount of oxygen in the system. 

In order to determine whether there are more stable 
ternary MS configurations we calculate a number of Li- 
Al-B phases with various stacking sequences and com- 
pare their relative stability, AEf, with respect to the 
phase separation in 11-LiB and O-MB2. The results in 
Fig. |4] suggest that 10-Li2AlB4 would be the most LiB- 
rich phase to form, while one could observe several AIB2- 
rich phases at different concentrations (a similar phase 
diagram is obtained for the Li-Ti-B system). The ef- 
fect of permutation of the "1" and "0" units on the en- 
ergy at given composition is illustrated with the 10- • -0, 
110- • -0, 11110- • -0 series shown in Fig. [U These re- 
sults clearly demonstrate the benefit of having stackings 
with alternating electron-poor "1"- and electron-rich "0"- 
units. For instance, 1010 is 20 meV/atom below 1100, 
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FIG. 3: (Color online). Top: location of known and predicted 
phases on the Li-Al-B phase diagram (the boron-rich Li-B 
and Al-B phases are not shown). Bottom: Relative stability 
of the proposed Li2AlB4 with respect to the decomposition 
along the lines shown on the top panel. 



which implies that the boron layer in the 11 block of 
1100, surrounded by two Li layers on each side, does not 
take full advantage of the available charge provided by 
Al. This behavior is consistent with our previous obser- 
vation that substitutional doping leads to the filling of 
boron states only in the layer closest to the dopant [24j. 

Note that the boundary of the convex hull MS2- 
LiB<-^C32-AlB2 is determined by the 10- • -0 series, which 
implies that there is a thermodynamic driving force 
for the formation of "1" -units in the "0" matrix. In- 
deed, doing extrapolation over the 10- • -0 series we find 
that replacement of one Al layer via rtAlB2 + 2Li = 
Li2Al„_iB2 rl + Al is an exothermic reaction with AE 
= 0.6 eV/(Li atom) (n — > oo). There is an alternative 
exothermic reaction with an even bigger AE — 0.8 eV/(Li 
atom) (n — > oo) that preserves the C32 structure: nAlB2 
+ Li = LiAl„_iB2 n + Al. Therefore, one would expect 
to see substitution of Al layers with Li first, and only at 
high Li concentration should the "l"-units start forming: 
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FIG. 4: Relative stability of Li2(i_a:) ALB2 ternary metal 
borides with respect to the phase separation into 11-LiB 
(MS2) and O-AIB2 (C32). The dotted lines connect members 
of the 110- • -0 and 11110- • -0 series. 



in the limiting case of 1:1 Li:Al composition the reaction 
LiAlB4 + Li = Li2AlB4 would be energetically favorable 
by AE — 0.4 eV/(Li atom). However, these reactions 
have apparently not been observed due to the formation 
of the stable ternary LiAlBi4 and the oxygen-containing 
phases 

In summary, based on our formation enthalpy calcu- 
lations, Li2AlB4 is marginally stable with respect to the 
considered known compounds under normal conditions 
but can be stabilized by hydrostatic pressure. Consid- 
ering the stabilizing effect of Al we expect a lower pres- 
sure threshold for the formation of the ternary MS-Li- 
Al-B phases compared to the case of the binary MS-Li-B 
phases. 

More detailed studies should be carried out in the fu- 
ture to account for the finite-temperature contributions 
in the Gibbs free energy; this will be a difficult task be- 
cause the boron-rich metal borides have large unit cells 
with fractional occupancies. In addition, one should con- 
sider not only compounds reported for a given ternary 
system, but also crystal structure phases observed in sim- 
ilar systems. For example, the absence of any stable Li- 
M-B phases (M = V, Nb, Ta) in the ICSD database 
14j could simply be an indication that these ternary sys- 
tems have not been fully explored experimentally yet. It 
would not be surprising then if attempts to synthesize 
the predicted layered phases would lead to the forma- 
tion of phases not considered here, such as the metal- 
rich Li2PcUB superconductor comprised of linear chains 
of boron [46(]. High-throughput simulation of selected 
ternary systems identified in our present work is subject 
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FIG. 5: (Color online). Band structure and partial density 
of states (PDOS) in 1010-Li 2 A1B 4 . The top panels show the 
region near the Fermi level (0 eV). 
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FIG. 6: (Color online). Band structure and partial density 
of states (PDOS) in 1010-Li 2 TiB 4 . The top panels show the 
region near the Fermi level (0 eV). 



of future studies. 



IV. Electronic properties 

It is interesting to see what changes in properties the 
combination of the non-superconducting MB 2 and po- 
tentially superconducting MS-LiB structures would in- 
duce. Electronic properties of the proposed MS ternary 
phases are examined for two representative 1010-Li2AlB4 
and 1010-Li 2 TiB4 compounds with hexagonal unit cells 
(1010 and 10 phases are nearly degenerate in energy). 
Boron layers in the 1010 structure are indistinguishable, 
but there are two different sets of degenerate boron states 
derived from those in the original "0" and "1" metal 
borides. For example, in MS2-LiB the practically dis- 
per signless va states of boron are 1.0 eV above the Fermi 
level 12|,l22j while in C32-A1B2 they are completely filled 



[.Ef =-1.6 eV, E£=-1.0 eV]; in 1010-Li 2 AlB 4 the two sets 
of boron pa states are not completely filled, being 0.6 
and 0.1 eV above the Fermi level. This is a very satis- 
fying outcome, because they still contribute a consider- 
able 0.042 states/ (eV-spin-boron atom) to the density of 
states (DOS) at the Fermi level (to be compared to 0.059 
and 0.049 states/ (eV-spin-boron atom) in MS2-LiB and 
C32-MgB 2 , respectively^]). Because boron orbitals do 
not overlap across the Li-filled portion of 1010-Li2AlB 4 , 
the pa states are not dispersed and give rise to the de- 
sired, nearly cylindrical two-dimensional Fermi surfaces 
similar to those of MS-LiB [H, El, In addition, 



the compound has pn Fermi surfaces, as the p7r-bands 



crossing now happens at about 0.7 eV. The p7r states 
DOS of 0.019 states/ (eV-spin- boron atom) is a signifi- 
cant improvement with respect to MS2-LiB, which lacks 
these states altogether. However, this contribution is still 
below the value in MgB 2 (0.032 states/(eV-spin- boron 
atom)) and, unfortunately, comes at the expense of los- 
ing some pa DOS. 

Figure [6] shows that boron in 1010-Li 2 TiB 4 is too over- 
doped to have any strong MgB 2 -type electron-phonon 
coupling. The bottom pa band at -1.5 eV is completely 
filled, and even though the top band is still 0.3 eV 
above the Fermi level the steep dE/dk derivative on 
the narrow cylinder along T-A results in a relatively 
low DOS of 0.010 states/(eV-spin-boronatom). 1010- 
Li 2 HfB 4 exhibits similar properties but has a slightly 
higher pa DOS of 0.017 states/(eV-spin-boronatom). 
1010-Li 2 MB 4 compounds with transition metals of va- 
lence V and higher present cases of borides with com- 
pletely filled boron pa states (see Table I) . 



V. Phonon softening 

As previously discussed [24| , the next step in the anal- 
ysis of superconducting properties of the proposed ma- 
terials should be an ah initio calculation of the electron- 
phonon coupling A. However, the full scale calculations 
become prohibitively expensive for compounds with large 
unit cells, and more so if they have nested Fermi surfaces. 
In the following we attempt to obtain information on the 
strength of the electron-phonon coupling from the com- 
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FIG. 7: (Color online). Frequency of the in- plane boron 
phonons at M (open symbols) and V (solid symbols) points as 
a function of boron-boron bond length. The known C32-MB2 
diborides are shown as blue squares, the proposed MS2-LiB 
is shown as red circles. The solid line shows a linear fit of 
the M-point frequencies for the known diborides with metal 
valence from III to V. 



FIG. 8: (Color online). Frequency of the in- plane boron 
phonons at M (open symbols) and F (solid symbols) for C32- 
MB2 diborides of group I, II, and III-V metals. The blue 
squares, black triangles, and red circles correspond to the 
known, unstable, and proposed (11-LiB and 1010-Li2-MB4) 
compounds, respectively. The linear fits to the sets of M 
point frequencies are labeled according to the valence of the 
metal. 



pounds' vibrational properties. 

Indeed, it has been previously argued that the unex- 
pected decrease of A in MS-LiB with respect to that 
in MgB2 can be related to the increase of the in-plane 
boron frequency and to the smaller softening of this 
mode. [23L \24 \ The considerations were based on the 



fact that for a zone-center optical phonon mode the 
softening due to screening by metallic electrons, Au>, is 
connected to the Fermi-surface average of the square of 
the electron-phonon matrix element, <?, through Au) = 
-4ui{g^}N(0)(N(0) is the electronic DOS at the Fermi 

level) ^SaEi- 

Unfortunately, there is no simple way to evaluate the 
unscreened frequency in real materials. However, in 
MgB2-type superconductors the M-point phonons have 
wavevectors larger than 2k p and are thus screened much 
less than the T-point phonons. For this reason we will use 
the ujm — wr difference to reveal how the in-plane boron 
vibrations soften across a set of binary and ternary lay- 
ered metal borides. It should be stressed that the above 
frequency difference cannot be treated as the real soften- 
ing for quantitative estimates of the electron-phonon cou- 
pling; nevertheless, we expect it to be useful for selecting 
candidates with the strongest electron-phonon coupling. 

For the C32-M B 2 diborides and MS2-LiB monoboride, 
we calculate the frequencies of the in-plane boron mode 



at r and the TO mode at M points in the frozen-phonon 
approach using the fourth-order corrections and the same 
settings as described in Refs. 12, [49|. We neglect the 
anharmonic and non-adiabatic effects since we are inter- 
ested in general trends. Full phonon spectra in existing 
metal diborides have been calculated previously 5^, 51 1 , 
and the general inverse dependence of the E^ g frequency 
at r on the boron-boron bond length has been pointed 
out [HH. In Fig. [7| we plot, as functions of the boron- 
boron bond length cib-b, the T- and M-point sets of 
calculated frequencies. The latter exhibits a much more 
consistent linear decrease with the bond length for all 
known metal diborides of groups III-V. One can further 
discern two subsets (Al, Hf, Zr, Y and the rest) which 
have T-point frequencies following the linear dependence 
with similar negative slopes (we see no apparent factor 
which could be responsible for this behavior). Note that 
for the electron-rich metal diborides ljm is below u>r, 
which illustrates why the former cannot be considered 
as the unscreened frequency. However, the calculated 
values of T- and M-point frequencies will be useful as a 
reference when we examine the vibrational properties of 
the ternary MS phases derived from the corresponding 
metal diborides. 

As expected, MgB 2 stands out as a compound with 
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a large softening of the frequency at T. In addition, 
Fig. [7] shows that the M-point frequency does not follow 
the general trend of group-III-V metal diborides either; 
this behavior cannot be directly related to the electron- 
phonon coupling with the pa boron states. The phe- 
nomenon is explained in Figure [5] where we plot calcu- 
lated frequencies for a large set of group-I-II metal di- 
borides. The M-point frequency ujm for divalent (Be, 
Mg, Ca) and monovalent (Li, Cu, Au, Ag, Na, K) metals 
shows similar and nearly perfect linear dependences on 
the boron-boron bond length. One can expect to have the 
same trends for the unscreened frequency, only it would 
be shifted up by a constant offset. 

These results give a new perspective on the relation 
between the vibration properties and the strength of 
the electron-phonon coupling in LiB discussed in Refs. 
[H, As a boride isovalent to MgB 2 , LiB follows 



the wm(ub-b) dependence obtained for the Be-Mg-Ca 
series of diborides shown in Fig. [H Therefore, one of 
the key factors in the LiB hardening of the in-plane fre- 
quency [with respect to that in MgB2] is the shortening 
of the boron-boron bond length, which depends on the 
metal and the particular filling of the pir and pa states 
of boron. Hence, it follows that the frequency softenings 
in LiB and MgB2 are, in fact, comparable. 

Finally, we can comment on the superconducting prop- 
erties of the proposed ternary compounds. The intro- 
duction of electron-rich metals (group IV and higher) 
into 1010-Li 2 MB4 overdopes the boron layers, and causes 
the vibrational properties to become closer to those of 
the corresponding diborides: ujt is still above u>m for 
Ti and Hf (Fig. |SJ). 1010-Li 2 AlB 4 is a more promis- 
ing MgB 2 -type superconducting material, judging by the 
pronounced softening of the in-plane boron mode (note 
that in A1B 2 o-r is harder than u>m)- Unfortunately, alu- 
minum induces the shortening of the bond length and 
a substantial increase of the in-plane boron frequency, 
which will likely weaken the electron-phonon coupling 
23L lil (the next section is devoted to a more accurate 
description of this compound). As for the proposed MS 
phases at other compositions, we find that all ternary 
systems with negative relative stability for 1010-Li 2 AfB4 
(Fig. [2|) have negative relative stability for 110-Li4AfBg 
as well (the latter are typically 2-7 meV/atom above the 
line connecting 11-LiB and 1010-Li 2 AfB4 and are only 
metastable). Therefore, if IIO-L14MB6 could still be syn- 
thesized with large in size transition metals, such as Hf 
or Ta, the superconducting properties of the boron layer 
in the "11" portion (which should be well isolated by the 
double layers of lithium from the electron-rich "0"-unit) 
would be enhanced due to the stretching of the over- 
all boron bond length in the compound. A similar phe- 
nomenon has been recently observed in stretched MgB 2 



films [52j. A recent Raman study of the vibrational and 
superconducting properties of doped MgB 2 has revealed 
that the adjustment of the fre quen cy alone might not be 
enough to achieve a higher T c [53[ . 

The trend shown in Figs. 0U is corroborated by 
the well known correlation between the increase of the 
electron-phonon coupling and the inevitable dynamical 
destabilization of the structure 54, 55, 56|, 57, 58| (a too 
strong renormalization can cause a Peierls-type distor- 
tion, a band Jahn- Teller transformation or a structural 
transition [55j). In our considered metal boride proto- 
types, the specific relation between stability and super- 
conductivity originates from the subtle effects associated 
with the filling of binding states of boron. 

The focus of our study has been the MgB 2 -type super- 
conductors, while there have been reports on Nb-deficient 
superconductors based on C32-NbB 2 with T c of nearly 10 
K [1, Q . Indeed, it would be interesting to explore the 
superconducting potential of Li 2 MB4 compounds (M = 
V,Nb,Ta) which would have a non-MgB 2 -type electron- 
phonon coupling mechanism. 



VI. Phonon spectrum and electron-phonon coupling. 

In the previous sections we have shown that Li 2 AlB4 
has a significant amount of boron pa and pit DOS at 
the Fermi level and exhibits a characteristic softening of 
the in-plane boron phonon mode at T. All these quanti- 
ties are lower compared to those in MgB 2 so that one 
can expect the superconducting properties of Li 2 AlB4 
to be less appealing. However, a full calculation of the 
electron-phonon coupling (for all the modes and all k- 
points in the Brillouin zone) is needed to understand the 
superconducting properties of Li 2 AlB4. For example, in 
MgB 2 the intercalant modes are weakly coupled; based 
on our results for LiB the situation could be differ- 
ent in Li 2 AlB4. In addition, due to symmetry constraints 
in calculations on periodic systems, energy minimization 
might not guarantee convergence to an equilibrium free 
of dynamical instabilities (imaginary phonon frequency 
at q 7^ r). Thus, dynamical stability has to be checked 
by explicit calculation of the phonon frequencies in the 
whole BZ. In this section we calculate vibrational phonon 
frequencies and the electron-phonon coupling of Li 2 AlB4 
using density functional theory in the linear response ap- 
proach [591 ] - 

We find that the theoretically devised structure of 
Li 2 AIB4 is dynamically stable with no imaginary phonon- 
frequencies. The obtained phonon density of states (PH- 
DOS) is plotted in Fig. [9l As it can be seen from the de- 
composition along selected cartesian vibrations, B modes 
are dominant at energies larger than 50 meV. In the low 
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energy region (< 50 meV), two clear peaks are seen, one 
at w 40 meV due to Li vibrations and the other one at 
30 meV due to Al vibrations. As expected from the large 
interlayer spacing between the B-layers, B-modes are flat 
in the direction perpendicular to the B-layers. 

The clcctron-phonon coupling A q „ for a phonon mode 
v with momentum q is: 

^ = 7, ATfn\ AT l5kri.k+qm| 2 ^( e kii)^( e k+qrn) 

(2) 

where the sum is over the Brillouin Zone. The matrix 
element is ff£„,k+qm = (kra|<5V/<5u q „|k + qm)/^2w q „, 



where 



is the amplitude of the displacement of 



the phonon, V is the Kohn-Sham potential and N(0) 
is the electronic density of states at the Fermi level. 
The calculated average electron-phonon coupling is A = 
Sqiy ^qu/Nq ~ 0.41 (Nk and N q are the k-space and 
q-space mesh dimensions, respectively [o0j|). 
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FIG. 9: (Color online). Phonon Density of States (PHDOS) 
decomposed over selected vibrations, Eliashberg function 
a 2 F(u>) and integrated Eliashberg function A(u>) in Li2AlP>4. 



The Eliashberg function 



1 



2N n ^ ' q 

qi/ 



(3) 



and the integral A(w) = 2 J Q W du'a 2 F(uj')/uj' are shown 
in Fig. [5] As can be seen, the contribution of the cou- 
pling due to in-plane vibration is substantially reduced 
with respect to MgB2. This confirm what could be qual- 
itatively inferred from the softening calculated in Fig. [8] 



On the contrary the coupling to intercalant modes is not 
negligible. 

The critical superconducting temperature is estimated 
using the McMillan formula [6l| : 



l c = exp 

1.2 P 



1.04(1 + A) 



A-//*(1 + 0.62A) 



(4) 



where is the screened Coulomb pseudopotential and 



/ \ v f + a 2 F(u) \os,(uj)/uj du 

{Lo)i og = e x Jo 



(5) 



the phonon frequencies logarithmic average. We obtain 
(ui)i g — 59.4 meV leading to T c of approximately 3.6 K 
for /i* = 0.1. This value might be enhanced by multiband 
effects. 



VII. Summary 

We have demonstrated that there is a number of Li-M- 
B systems, in which MS ternary borides gain in stability 
with respect to MS2-L1B and C32-MB 2 . The finding 
suggests that these potential superconductors could be 
synthesized under normal conditions or grown on metal 
diboride substrates with a matching lattice constant. We 
show that some boron states in the MS ternary com- 
pounds [derived from those in MS2-LiB and C32-MB2] 
are still hole-doped and could give rise to the MgB2-type 
superconductivity. In order to pre-select compounds with 
the strongest electron-phonon coupling we examine the 
softening of the in-plane boron phonon mode in a large 
class of metal borides. We find a very well defined corre- 
lation between the frequency of this mode and both the 
boron-boron bond length and the valence of the metal. 
This observation allows us to give a new interpretation of 
the hardening of the in-plane boron phonon mode in MS- 
LiB and identify one compound, Li2AlB4, that should be 
a superconductor with T c of about 3.6 K. 
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